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A contribution to 7/\|/ hadroproduction is analyzed in which the meson production is 
mediated by three-gluon partonic state, with two gluons coming from the target and one 
gluon from the projectile. This mechanism involves double gluon density in one of the pro¬ 
tons, hence this conhibution enters at a non-leading twist. It is, however, relevant due to 
an enhancement factor coming from large double gluon density at small x. We calculate 
the three-gluon contribution to //\|/ hadroproduction within perturbative QCD in the kr- 
factorization framework. Results are obtained for differential pT'-dependent cross-sections 
for all 7/\|/ polarizations and for the sum over the polarization components. The rescatter¬ 
ing contribution is found to provide a significant correction to the standard leading twist 
cross-section at the energies of the Tevahon or the LHC at moderate pj. We suggest J/^^f 
production in proton-nucleus collision as a possible probe of the triple gluon mechanism. 


I. INTRODUCTION 


Production of heavy vector quarkonia in proton collisions has been a subject of intense experimental 
and theoretical investigations since CDF found huge excess of the measured cross-sections of prompt 
ttlD with respect to the standard theoretical predictions in QCD. Initially the theory predictions 
were based on a leading order collinear QCD approximation, in which the hard matrix element of a 
partonic color singlet subprocess gg —?■ J/'^g gives the dominant contribution. This approximation, 
however, underestimates badly the measured cross-section and yields an incorrect pp distribution of the 
produced mesons, so it was necessary to consider alternative mechanisms of quarkonia hadroproduction. 


In p^icular, two theoretical concepts received a lot of attention: the so ca’ 
and the other approach based on the ky-factorization (KTF) scheme 


led CO 


or octet model (COM) 


IttilSll. For a comprehensive 

theory review of various approaches to quarkonium hadroproduction see o. 

The COM is based on the fact that the hadron wave functions contain higher Fock components. For 
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heavy quarkonia it implies that the heavy quark-antiquark (QQ) pair in a heavy meson may be aeeom- 
panied e.g. by one or more gluons. So the QQ-pair may appear in the meson in a eolor oetet state yfl. 
Amplitudes of such octet components are subleading (w.r.t. the standard, color singlet component) in 
heavy quark velocity expansion within the Heavy Quark Effective Theory OO, but in the vector meson 
hadroproduction the color octet states contributions are enhanced by the corresponding hard partonic 
cross-section for production of such states oildl]. With only a few free parameters the color octet model 
is quite successful in describing the magnitude and pr-dependence of the prompt 7/\|/ hadroproduction 
[i however, it faces some problems with an accurate desription of the produced mesons polarization. 
In p^icular, predictions for the meson polarization change a lot when going from LO to NLO accuracy 


fl 


90. Besides, neither LO nor NLO COM estimates of prompt J/xf/ polarization are fully consistent 
with the data. It is fair to add that in the COM the amplitudes of color octet QQ transitions into mesons 
are not known from the first principles, and the theoretical estimates of the cross-section magnitudes rely 
upon fitted parameters. On the other hand, kinematical dependence of the cross-sections, e.g. the pr- 
dependence, is already a genuine QCD prediction (for each partonic subprocess separately) and the good 
description of the meson pr-dependence within the COM supports strongly validity of this approach. 
The COM may be also understood within hard factorization theorem in QCD, where universal (inde¬ 


pendent on the environment) fragmentation functions exist o 


numbers (including color octet) into final state quarkonia [Il7ll . 


QQ partonic states with various quantum 


Another theoretical approach to quarkonia hadroproduction relies on the ky-factorization scheme. 
In this framework, one allows for non-zero transverse momentum of colliding partons. Because of 
multiple emissions in the course of QCD evolution the parton momentum may become sizeable if the 
evolution length is sufficiently large. This effect is especially pronounced if gluons are probed at small x, 
which is the case for prompt quarkonia production at the Tevatron or at the LHC. Then the kinemati¬ 
cal distributions of produced quarkonia may be strongly affected by kr of incoming gluons. Indeed, 
inclusion of gluon kj in the analysis brings theoretical estimates much closer to the experimental data, 
both for the magnitude and for the shape of the distributions, already at the lowest order in QCD. This 
works rather well for the C-even charmonia, like Xc, but for the accurate description of Tevatron data 
on p rompt 7/\|/ and \|/' production a color octet component was necessary also in the KTL approach 
lll2ll . although, smaller than in the collinear approximation. Lor the recent LHC data on prompt char¬ 
monia production, however, a good theoretical description of magnitude, kinematics and polarization 
was found within KTL without including the color octet mechanism [Il4ll . This good description in KTL 
CSM was achieved by using a more advanced unintegrated parton density, CCLMAO lllSn . Still, it is fair 
to say that a global, consistent description of quarkonia production in KTL has not been formulated yet. 
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It follows from this discussion that theoretical understanding of heavy quarkonia hadroproduetion is 
not yet fully satisfaetory. Polarization deseription and the need of fitting key parameters are somewhat 
weak points of the eolor oetet model, and a global deseription of 7/\|/, \|/' and %c hardroproduetion was 
not formulated yet in the KTF approach. Therefore it is possible, that both proposed approaches are not 
eomplete and need to be supplemented by another meehanism. 

The COM and CSM approaehes both in the eollinear and KTF framework share an important eom- 
mon feature — it is assumed that at the parton level the quarkonia produetion is initiated by two gluons. 
At very large energies, however, when the gluon density beeomes large, it is expeeted that sizable eor- 
reetions to this pieture may eome from subprocess with more than two initial state partons, which would 
be elassified as reseattering or shadowing/anti-shadowing eontributions. Sueh effeets are expeeted to be 
partieularly important in eollisions with heavier nuelei. First estimates of sueh multiple-gluon eontribu- 


tions to heavy c 


uarkonia produetion in pp and pp eollisions were performed by Khoze, Martin, Ryskin 


and Stirling [Il9n . who eonsidered three-gluon initial partonie state. Reeently, Ma and Venugopolan [|20n 


prop osed resummation of reseattering eontributions within the eolor glass eondensate (CGC) formalism 
2311 . A good deseription of unpolarized data for y/\|/ and \\f' was aehieved down to very low pr by 


m 


matehing the resummed CGC eross-seetion with predietions of the eollinear NRQCD COM approaeh 


m 


2M- 


In this paper, following the idea of Khoze, Martin, Ryskin and Stirling [Il9h . we shall investigate the 


three gluon meehanism of quarkonia hadroproduetion in whieh hard reseattering plays the key role in 
eolor neutralization of the QQ pair. Thus, we shall eonsider the 3-gluon fusion proeess at the parton 
level: 3g —y/\|/ in the eolor singlet ehannel instead of the standard eolor singlet proeess gg —?■ J/'^g. 
This three gluon fusion partonie proeess is, of eourse, a higher twist eorreetion to the eross-seetion, as 
it eouples to a double gluon density in one of the protons. Therefore, one expeets it to be suppressed 
w.r.t. the standard, leading twist eontribution by powers of the hard seale. In high energy eollisions, 
however, this reseattering term reeeives signifieant enhaneement due to large gluon density, refieeted 
by large ratio of double and single gluon density at small-x. Moreover, the matrix elements of the 
reseattering term 3g J/\\r, and the standard one gg —> //t|/g, eome at the same order of perturbation 
theory, so no additional suppression of the reseattering pieee oceurs in the matrix element. Henee, 
despite of the higher twist nature, it may be still important. As already stated, sueh proeesses were 
already proposed in Ref. [Il9ll where estimates were given showing that reseattering might explain the 
diserepaney between the data and the predietions based on the eolor singlet meehanism. We address and 
develop this idea in more detail, by performing a eomplete ealeulation of the reseattering eontribution 
in the KTF framework, ineluding an explieit ealeulation of pr distribution of produced quarkonia and 
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FIG. 1: a) The amplitude M of the 7/\|/ hadroproduction in the triple gluon fusion mechanism; b) the dominant 
diagram in the square of the amplitude, \9v[\'^. 

the polarization composition. We find that the rescattering correction is a significant contribution to the 
prompt J/'^ cross-section especially at moderate pr and should be included in a precise description 
of the data. It should be stressed that the rescattering correction complements in natural way both the 
COM and the standard KTF approach, providing a kinematics- and process-dependent mechanism of 
color neutralization. Also, we expect that the rescattering mechanism should be strongly enhanced and 
thus - even more important for heavy quarkonia production in collisions with nuclei. 

The paper is organized as follows: in Sec. HI] we introduce the ki nematics of the process and the 
notation, in Sec. imi we derive analytic formulae for the amplitudes and cross-sections of the direct 
quarkonium production by rescattering, in Sec. |IV|we give numerical estimates for the cross sections 
and conclusions are given in Sec. jV] Technical pieces of derivations are given in the Appendix. 


II. KINEMATICS 

We consider direct inclusive production of heavy quarkonia in high energy proton-(anti-) proton 
collisions and focus on a process mediated by partonic subprocess with tree incoming gluons: 3g —)■ V. 
In general, V can be any heavy vector quarkonium, \\f', T and so on, but we choose the best 
measured J/\\f production as the reference process. The three-gluon fusion contribution enters as a 
higher twist correction to the direct vector quarkonium production. The notation for the four-momenta 
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of incoming protons {pA, Pb), the outgoing meson, p, outgoing hadronie final states for A' and p'g 
for B') is explained in the left panel of Fig. [TJ 

We analyze the process within perturbative QCD in the high energy approximation, using the frame¬ 
work of kf-factorization. In the high energy limit, only the leading power of 5 = {pA is retained. 

In this approximation the proton masses may be negleeted, p\= Pb = 0, and s ^ 2pA- Pb- The masses 
of outgoing states are denoted by p^ = p^, = M^,, = M^, for 7/\|/ and hadronie final states A' 

and B' emerging from the target and the projeetile eorrespondingly. 

The diagram in Fig. [T^ illustrates the partonic topology of the amplitude for direet vector quarko- 
nium production with three intermediate gluons. The four-momenta of intermediate gluon eouplings 
to Pa —> a' and pb ^ B' transitions are denoted by k\, ki and I eorrespondingly. The four-momentum 
eonservation imposes eonstraints on the four-momenta, 

PA=PA' + h+k2, Pb=Pb' + 1, P = ki+k2 + l. (1) 


We employ the Sudakov parameterization of the four-momenta. 


PA'{B') = <^A{B)PA + ^A{B)PB + Pa'{B')L , P = apA + ^PB + Pi ., 


( 2 ) 


where the transverse direetions are in the plane orthogonal to the direetion of the eollision axis in the 
CMS frame. Then, the kinematie eonstraints lead to: 


2j_n/r2 „2 I 11^2 I «/f2 


PA = 




= 


, |3 = 


(3) 


The proeess of 7/\|/ produetion is analyzed in the kinematieal domain where = 0(1), (3 b = 0(1) and 
consequently (3^ ~ ~ 0(1/5). Thus, the gluon momenta satisfy 


k = ki+k2 ^ {l-aA)pA-PA'±, I ~ {^-%)pb-Pb'. 


(4) 


and from the last term of O one gets 


a^il-aA, (3?«l-ps, p^^ -pa'a-Pb'l- (5) 

In the Sudakov parameterization of the gluon momenta k,- = a/pA + 13fps + k,j^, I = d^pA + (3^Pb + 
(see Fig.d^) one has relations: 


a^-fa 2 = a, + = 0{l/s), d = 0{l/s), (3i = (3. 


( 6 ) 
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III. THE TRIPLE GLUON CONTRIBUTION: ANALYTIC FORMULAE 

A. The triple gluon amplitude 


The amplitude in Fig. [T^ reads 


1 h 


4-^aiplp2 ’ ^2, F)- ^2l2 - -^Aviva fci , k2) 


(271) 


/C 1 /C 2 


(7) 


where = (pfi',/|7ji'(0)|pz„rB) and ^ 1 ,^ 2 ) give the amplitudes of finding gluon with 

eolor b 1 , four-momentum I and two gluons of colors a 1 , a 2 and momenta , ^2 = ^ — respectively in 
the decay products of a struck proton. The incoming proton polarizations are denoted by rA,s. The vertex 
SoJpJp 2 (/,ki,k 2 ;p,e) describes the amplitude of 7/\jt production with momentum p and polarization 
£. In the standard Regge kinematics of high energy scattering at a small momentum transfer one can 
approximate polarization tensors of gluons (in the Feynman gauge) as and ~ 

—2p^p'g/s for i — 1,2. Then, one can write the amplitude (Q in the following form 




ki{k — kiY 




( 8 ) 


where the impact factors (at the amplitude level) for a pg —5' -f g, pa —^ + 2g and 2g -|- g —)■ 7/\|/ 

transitions read correspondingly: 


^‘'j^''^{pA,rAA'\^k,^k'M,k-kx) = J d^\s‘^^^^^^{pA,rA,A' ,h,k - ki) 


p'sPf 


Ol „Pl „P2 




(9) 

( 10 ) 

( 11 ) 


The total contribution of three gluons to 7/\|/production amplitude comes from a diagram depicted in 
Fig.[Bi and an ‘upside-down’ diagram with the same topology, but with two gluons coupling to the lower 
vertex pb ^ B'. The contribution of the other diagram, M' may be obtained from OvC, corresponding 
to the diagram in Fig. [T^, by exchange of kinematical variables: a -fA |3. It is important to note that the 
interference of these two amplitudes, M and M' is a subleading effect. This is because terms fM*!M 
and are driven by the QCD evolution of a four-gluon 7-channel state (two ladders in the large 

Nc limit) emerging from one proton and a two-gluon 7-channel state from the other proton (see Fig. 
[2^) whereas the interference terms and are driven by subleading three-gluon 7-channel 

states from both sides, see Fig. 125. Hence, at large energies and for large scales, the interference terms. 
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FIG. 2: Topologies contributing to the rescattering correction in the 7/\|/ cross-section : a) the topology with three 
two-gluon ladders (leading), b) the topology of two three-gluon BKP states (suppressed). 

shown in Fig.[2j3, may be neglected and the contribution to the J/'^ production cross-section from three 
intermediate gluons do \ 


B. The uncorrelated triple gluon cross-section 


The simplest model of the two-gluon distribution in the proton assumes lack of correlations of the 
gluons in the transverse plane. This means that the double gluon density is proportional to a prod¬ 
uct of independent single gluon distributions. Below we estimate the triple gluon contribution to 7/\|/ 
hadroproduction in this scenario. 

The contribution to 7/\|/ cross-section, corresponding to the diagram shown in Fig. [IJ 5 takes the 
following form: 


do 




/ 


dH 1 


X 


{2nfJ (27t)4(/2)2 2 ,p 
d^k r d^ki f d^k 




( 12 ) 


f d k f d k\ f 

J J kf{k — ki)^J 


:4>' 




(ttfe, pp, A:i, A: - A:i; A:3, A: - A:3) 


X 


,k-ki, 

d^p 


kl{k-k,Y 


X (27tr5(Z-fA-p)^0(po)27t5(p^-Mj/^) + (t/oL), 


where the last term describes contribution from the graph where upper and lower vertices from Fig. [IJ 5 
are interchanged. The pomeron - proton impact factor is defined by the formula 


2,p 


(a/,p/,/) = ^ / J 4 >p/( 27 t)'^ 5 (/-ps + ps/)S 

^ rs B' •’ 


b\ 

Spi 


^b2* Pa Pa 


(13) 
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ei,6i,/3i,/i ^ e2,h2,^2-,h 



FIG. 3: The unintegrated gluon distribution. 


whereas 

= !d^A'{'^T^Y^{k-pA+ pa') (14) 

^ rA A' 

X ^‘2‘'^{pA,rA,A'{pA.rAA'-k^) 

gives the two-pomeron-proton impact factor (2P-p). In the last line the energy-momentum conservation 
delta function in the 7/t|/ vertex gives a, (3; (3. 

Several remarks are in order here: 

1. The pomeron - proton impact factor d> 2 ,p (see Fig. [3]) integrated over longitudinal variables gives 

j da,d^„^\'^(a,Ai,l) = (15) 

where P/ ~ |3 follows from the energy-momentum conservation in 7/\|/ vertex and /(P,/^) is an 
unintegrated gluon distribution. 


2. The 2P-p impact factor integrated over longitudinal variable P^t can be decompose in the following 
way (see Appendix for the derivation): 

I dakdhK,r'"'^^k.hM.k-ki-,k3.k-k3) = 

p2 (2'Ji)'^ 

,.g* , 2 -*3) 

(^c - 1) ^ 

+ 5“'“45fl2«3y(oc,k?)/(a,(A:-A:i)2)S^(A:i+A:3-A:) + 5"'"25«3«4___] ^ (i6) 


where ~ a follows from the energy-momentum conservation in y/\|/ vertex. The last term in 
parenthesis does not contribute to the cross-section due to the color factor contraction with the 
y/\jt vertex. S(k) is a symmetric function peaked around A: = 0 and Rgh is the Shuvaev factor I125n 
introduced to account for the fact that longitudinal variables of the gluons are not equal and the 
off-diagonal gluon distributions enter the cross-section. 
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3. The y/\|/ triple gluon vertex reads ll26ll : 




JN 


^aya2b\ 

Nr 




V>/Y(a,(3;/,A:i,A:2;e) = 47tmcgy/x,, 


£*-{xpA + l±) 


+ 


£ -PaKX- 


r+ {ki-^k2Y+ Aml 

Aki-k2' 


ys 


’-)+£* ■Ia_ 


N + {ki -A:2)^ + 4m2 


(17) 


where k 2 =k — ki. 


Substituting equations (fT3l) through (fTTI) into (fT2l) performing delta funetion integrations and eolor 
faetor eontraetions one arrives at the final formula for the reseattering eorreetion to polarization depen¬ 
dent differential J/\\r hadroproduetion eross-seetion, 

^ 2 N^-4 Rlt, ct;/(P, (p~kf)f(a,k])f(a,[k-kif) 

X |V}/^(a,|3;A:i,A:-A:i,/7-A:;£)|^-h(a^|3,pA^Ps), (18) 

where the last term deseribes eontribution from the graph as in Fig. [TJ? where four gluons are attached 
to the pb instead of pA vertex. In the derivation one needs to evaluate an integral f d^kj_S^(kj_) which, 
assuming the Gaussian form of S(b), may be related to the inverse of multiple scattering parameter, 
fd^k^S^(k^) = 


C. Triple gluon cross-section with gluon correlations 


Formula (fTSl) describes the distributions of transverse positions of two-gluons coming from one of 
the protons as being independent, i.e. uncorrelated, in the transverse plane. At large pj of the meson 
one should, however, expect some correlations of these distribution to emerge in the course of QCD 


evolution of double gluon density. In this evolution 11271] . a single parton ladder may split into two 


ladders, emerging at the same impact parameter, that introduces correlations of the parton positions. 

A diagram that describes the gluon parton splitting into two gluons is shown in Fig. 01 This contribu¬ 
tion to hadroproduetion was discussed in lll9l] in the collinear limit and using only an approximate 
estimate of QCD evolution of the double gluon state. This correlated reseattering correction is po¬ 
tentially important as it is expected to lead to a less steep pr dependence of the resulting differential 
cross-section component than the one obtained with the uncorrelated gluon distributions. We calculated 
this amplitude with full kj dependence of the three gluons that enter the J/\\f vertex, an assuming that 
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FIG. 4: The correlated triple gluon contribution. 

the initial gluon in a upper vertex, that acts as a source of the double gluon, is collinear with the parent 
proton. After standard steps we obtained the following cross-section, 

A/^-4 r r cf-kcfkicf-k 2 f{^,(p-kf-) 

d\n^dp\ AN^{Nc—\-)J^ J {p—kYk\{k — k\Yk2{k — k2)'^ 

X Vji^{a,^\ki,k-ki,p -k-,z)V]i^{a,^\k2,k-k2,p -k\z). (19) 

An analogous formula is valid also for a quark from the target replacing the gluon as a source of the two 
correlated t-channel gluon after the suitable replacement of the pdfs and the color factor. This formula 
provides the lowest order approximation to the correlated rescattering amplitude, which is only an input 
amplitude for the evolution of the two-gluon state from the parent gluon to the y/\|/ vertex. Since the 
typical evolution length in rapidity is sizable for this evolution one expects a potentially large effects of 
the QCD evolution to occur. The treatment of this evolution in the framework of the small x resummation 
is possible but it is a highly nontrivial task and it is beyond the scope of the present analysis. However, 
we estimated numerically the lowest order contribution of (fT^ and found it to be significantly smaller 
than the uncorrelated rescattering contribution of (fTSl) . Thus in the present analysis this contribution will 
be disregarded. 


IV. NUMERICAL RESULTS 


A. Triple gluon corrections in pp and pp collisions. 


In the numerical eva 


uations the unintegrated gluon densities were used derived from the CTIO 


collinear gluon density OOn using Kimber-Martin-Ryskin (KMR) scheme 113 Ih with the factorization 
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a) b) 


FIG. 5: Differential cross-sections for 7/\|/ hadroproduction: a) 


dpT 


iron, = 1.96 TeV and b) 


da 


|y|<o.6 


X Br(7/\|/- 




at the Teva- 




•d¥-dpT 


X Br(7/\|/ —)• u+n ) at the LHC, Jl = 7 TeV. The data-points shown: 

lr|<0.75 

a) CDF Il28ll and b) ATLAS Il29ll . The curves shown: bold continuous the triple gluon contribution computed 


in this paper, the dashed and dotted thin curves represent results of Ref. 


in the collinear approach: color singlet 


(LO and NLO) and combined color octet and color singlet results (see the legend in the plots). 


scale given by the transverse meson mass, ^ + Pj- The running strong eoupling eonstant of a 

gluon with virtuality was evaluated at the seale + k^, with Me — This seale choice 

is eonsistent with the KMR preseription. We set the multiple seattering parameter value Oeff =15 mb, in 
aeeordanee with the experimental results from the Tevatron 11321 1331] . For the effective gluon distribution 
X in the Shuvaev faetor R^h we set X = 0.3, whieh leads to Rsh = 1-3. 

In order to visualise the relevanee of the triple gluon correetion we eompare the obtained results to 
a limited ehoice of experimental results from the Tevatron and the LHC, and to the standard collinear 
QCD fits within the eolor singlet and oetet models. In this paper we do not aim to fit the data, nor provide 
the global description, so for clarity, only px dependencies near the central rapidity 7 = 0 are shown. 

In Fig. [5^ the results of numerieal evaluation of the triple gluon eorreetion to the differential meson 


produetion eross seetion 


da 


'pp^J/^VX 


dpT 


|y|<o.6 


X Br(y/\|/ ^ ) (ineluding the branehing ratio of the me¬ 


son decay to muons) are shown for the Tevatron energy (^/s = 1.96 GeV) — the bold continuous line. 
For reference we display also the CDF data 11281] and the results of eollinear ealeulations of the eolor 
singlet and CS-i-CO predietions at LO and NLO [a]. Clearly, the triple gluon eontribution enters as a 
subleading eorreetion to the prompt y/\|/ eross-seetion that may reaeh 20-25% of the total eross-seetion 
at moderate pp and beeomes negligible at larger pp. On the other hand, the triple gluon eontribution 
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exceeds the standard (collinear) color singlet contribution in the relevant pj range. 

In Fig.[5jr we also show a similar comparison for the ATLAS data in the central rapidity region. Note 

X Br(7/\|/ 


da 




|K|<0.75 


that the data-points describe here the double-differential cross-section, 

The overall pattern of different contributions and the data-points is similar to the pattern de¬ 
scribed above for the Tevatron energies. Again, the triple gluon contribution exceeds the CSM contri¬ 
butions, but makes not more than 20-25% of the total prompt 7/\|r cross-section, and the triple gluon 
correction becomes negligible at larger pj. 

In Fig. polarized components are shown of the triple gluon correction in the helicity frame for the 
Tevatron (at y/s = 1.96 TeV) and in Fig. [ 6 b for the LHC (at 1/5 = 7 TeV) and for the central rapidity, 
y = 0. The transverse polarization 1 in the plot is referred to the transverse polarization contained 
in the plane spanned by the beam axis and the meson three-momentum in the laboratory frame, and 
the transverse polarization 2 is perpendicular to this plane. Clearly, at low transverse momentum the 
longitudinal and the total transverse cross-sections are close to each other, and with increasing pj the 
longitudinal component becomes dominant, almost saturating the total triple gluon correction at px = 20 
GeV. As seen from Fig. the pattern of polarized cross-section for the Tevatron and the LHC is very 
similar. 

The presented results for the triple gluon correction in 7/\|/ hadroproduction are our central theoretical 
predictions. There is, however, theoretical uncertainty of these results coming from scale choices in the 
running coupling constant, the details of the unintegrated gluon densities, and from the unknown higher 
order corrections. The full analysis of these uncertainties is beyond the scope of this paper, however, 
we performed some first estimates and tested sensitivity of the results to the pQ scale variation in (1s{pq) 
and collinear parton choice in the KMR scheme. The results depend weakly on the collinear parton 
set, however, the increase of pQ to IpQ leads to reduction of the cross-sections by about 50% in the 
moderate pj range. For increasing pj this effect is getting weaker. On the other hand, applications 
of the unintegrated gluon densities (and related schemes) CCFMAO lilSi] and JH2013 ll34l] instead of 
the KMR scheme gave results larger by a factor of 2-3 in moderate p^-s (say for pr <5 GeV), and at 
larger py-s (say for pr > 7 GeV) the CCFMAO and JH2013 results are close to our central predictions. 
Summarizing these theoretical uncertainty checks, the normalization of our predictions in uncertain by 
a factor 2 up or down for moderate pr-s, and the normalization and shape of the distributions at large 
pr is much less uncertain. 





13 




a) 


b) 


c) 


FIG. 6: Polarization components of the triple gluon correction: a) (f'Opp^jdPT\y_Q for the Tevatron 
at ^/s = 1.96 TeV, b) LHC at y/s = 7 TeV, c) ratio of the cross-section for 

longitudinal and transverse polarizations. In a) and b) the polarizations a are chosen in the helicity frame. 

B. Probing the triple gluon contribution with nuclear beams 


In pp and pp the contributions of triple (or multiple) gluon exchange mechanisms of heavy vector 
meson hadroproduction may be not easy to disentangle from the standard mechanism of production 
within the CSM and COM. In particular, the triple gluon contribution may be to some extent absorbed 
into fit parameters of the COM. A more sensitive probe of the triple and multiple gluon exchange effects 
in vector quarkonia hadroproduction should be provided in experiments with nuclear beams, in particu¬ 
lar, in proton-nucleus collisions, pA, where A is the number of nucleons in the nucleus. In such collisions 
the triple gluon contribution scales with A in different way than the standard (two-gluon) contributions. 

Let us write the pA cross-section as a sum of the standard gluon-gluon contribution Ogf, a 
contribution with two gluons coming from the proton and a contribution with two gluon coming 


from a nucleus o 


pA 

sigsV 








( 20 ) 

In particular, for the pp scattering it reduces to OyP = Ogg + ^‘{^g-^g + ^^g[gg) ■ Since in Ogg and in 
one probes a single density of gluons in the nucleus, these cross-sections are enhanced by a factor of A 
w.r.t. the proton-proton cross-sections, Ogg =AOgg and 

The triple gluon contribution to 7/\|/ production with two gluons coming from the nucleus is en¬ 
hanced by A^R^/R\, with Rp and Ra being the proton and the nucleus effective radius. This dependence 
on the radius Ra comes from the integral of the nuclear impact parameter density Sa(6 ) of the gluon 
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distribution, / d^bS\{b) ~ ^/R\- Since Ra ^ A^/^Rp, we have Combining the 

A-dependencies of the chosen cross-section components one finds for the nuclear modification factor, 

^PP 


RpA = 


pA 

Aof 


= 1 + (A^/3-1). 


sigg) 

^pp 


( 21 ) 


Thus the deviation of RpA from one, 5^ = RpA — 1, measures the nuclear effects of multiple gluons from 
the nucleus. In our approach we neglect the effects of more than two gluons coming from the nucleus 
which is only motivated if the triple gluon correction is not large, 5 a 1 . 

Based on the obtained results on triple gluon contribution we also performed an estimate of 5 a for 
ALICE experiment conditions, that is for y/tjt inclusive production in pVh collisions at y/s = 5.02 TeV, 
-4.46 <Y < —2.96 (the nucleus fragmentation region) and 2.03 < 7 < 3.53 (the forward region). We 
found 5 a — 20% for in the backward region and 5 a — 100% in the forward region. The result for 
the backward region is slightly larger from the ALICE measurements showing evidence for about 10% 
nuclear enhancement of RpA with experimental errors of about 10% 11351] . Thus our estimates of the triple 
gluon correction are at the upper limit of ALICE measurements and more precise measurements may 
be used to constrain the details of our calculation and of the unintegrated gluon densities. The available 
data prefer a conservative approach of our central theory prediction with the unintegrated gluon densities 
derived from the collinear gluon densities within the KMR scheme. 

In the forward region the triple gluon correction to RpA is found to be large, over 100%. This implies 
that higher order rescattering corrections, with more than two gluons from the nucleus, are also important 
and the triple gluon contribution does not provide a reliable estimate of nuclear effects in the 

forward region for pPb collisions at the LHC. Therefore, in this region one should perform resummation 
of multiple scattering effects, as it was done e.g. in the color glass condensate approach |120il2^. We 
conclude that the triple gluon contributions to 7/\|/ hadroproduction may be experimentally constrained 
in pA measurements with lighter nuclei and/or central to backward meson rapidities where the triple 
gluon contribution is expected to bring a dominant correction to RpA- The key observable for such 
measurement is the A-dependence of the nuclear modification factor RpA- 


V, CONCLUSIONS 

The performed calculations show that the uncorrelated triple gluon contribution to 7/\|/ hadropro¬ 
duction introduce a sizable ~ 20% correction to theoretical predictions of prompt J/\\r hadroproduction. 
The relative correction is largest for moderate transverse momenta pr < 8 GeV and it quickly becomes 
negligible with pp growing beyond this range. The result for moderate pj is uncertain by about factor 
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of 2 (up or down), but the large pj shape and normalization is more stable against variations of the model 
details. It follows, that the uneorrelated eolor singlet triple gluon eontribution is unable to explain the 
large exeess of the 7/t|/ hadroproduction eross-section over the collinear eolor singlet predictions, espe¬ 
cially at larger values of pr- The triple gluon contribution, however, is found to be larger than the NLO 
color singlet cross-sections and therefore it should be relevant for theoretical analyses of heavy vector 
quarkonia hadroproduction. In particular, the estimated triple gluon contribution exhibits strong pj de¬ 
pendence of the meson polarization. At small pr this mechanism leads to equal rates of longitudinally 
and transversally polarized mesons, whereas at larger pj the longitudinal polarization strongly domi¬ 
nates. So the rescattering correction may modify the polarization pattern following from the collinear 
COM or KTF descriptions and it should be advantageous to complement with this contribution the fits 
to quarkonia hadroproduction cross-sections. 

For proper theoretical understanding of heavy vector meson hadroproduction it should be important to 
constrain experimentally the triple gluon (or rescattering) corrections. This may be not easy in pp and pp 
collisions where free parameters of the collinear COM or KTF fits may be used to absorb this correction. 
The multiple gluon effects, however, are characterized by a distinct (nonlinear) dependence of the gluon 
densities of colliding beams. Such effects may be directly probed in proton-nucleus {pA) collisions, 
where the relative triple gluon correction should be enhanced by a factor of A^/^. Recent ALICE data on 
inclusive 7/t|/ production in pPb collisions at i/s' = 5 TeV in the nucleus fragmentation region show some 
evidence of such non-linear nuclear enhancement, however, more extended and accurate measurements 
are needed to determine experimentally the magnitude of the triple gluon correction with satisfactory 
precision. 
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Appendix A 

The relation between 2P-p impact factor and unintegrated gluon distributions is read out from the 
collinear limit of semi-inclusive two particle production process pp ^ V 1 V 2 X described by the amplitude 
given in Eig. |7J The total cross-section in the impact parameter space is given by the collinear formula 


16 



FIG. 7: The amplitude for a semi-inclusive two particle production process pp —)• V 1 V 2 X 


(Fig. lit) 


a = 6i(ai,pi)62(a2,p2)g(oci,A^)g(oc2,A^)g(Pi,A^)g(p2,A^) 

X jd^bd^b^dhis{bi)S{b-bi)S{bi)S{b-bi), (Al) 

where we assumed that 2 -gluons density distribution faetorizes into two single gluon distributions g(a) 
and 61,2 are eross-seetions for the gg —)■ Vi ^2 produetions at the partonie level. The funetion S{b), whieh 
is approximated by the Gaussian form, S{b) = e'Kp{—b^/R^)/nR^, is normalized 

J d^bSib) = 1, 

where Rp is the effeetive proton radius. In the momentum space the total cross-section reads 

Oeff 

where 

and the normalization condition translates into S{0) = 1. 

The partonie amplitude for the production gg V\ has a form (see Fig. [ 8 b) 


where describes a production amplitude and is a gluon polarization vector. The cross-section, 

in leading 5 approximation, reads: 


di = 


715 aiPi 

- 1 k\J\ 


|Fi(fci,/i,pi)|^5(ai|3i5 



(A4) 
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FIG. 8: Colliding hadrons in the impact parameter space (left). The partonic amplitude for the production process 
gg V'l (right). 


where k\ a\pA+ki^,li ^ ^[pb + h±,Pi = clipa + ^iPb + Pi± and ai ttj, Pi Pj. The effective 
vertex is related to the production amplitude via 

The amplitude for the process in Fig.|7]reads 


-iM = 


J J^^Av%iPA.rAA\kuk2) ^lp^ai(^FPl -kuPl) 




and k 2 = Pa— Pa' — k\■ The cross section is then expressed 


as 


- 1 

dOpp^ViV2X = J7- = L 1 ^ 1 ^’ 

HPA-Pb) ^rA,rBA'.B' 

and the phase space integral takes the form 


1 


(A6) 


d^ 


‘pp^ViViX 


= {2Tlfh{pA+PB-PA'-PB'-P\ - P2)d^A'd^B'd^pid^P2- 


(A7) 


The integrals d^A'-id^p' describe the phase space integrations of A' and B' final states whereas 
, d^p^ correspond to single-particle phases-spaces for the Vi, V 2 particles. 

Using the appropriate eikonal approximation for gluon polarization tensors and inserting unities in 
the form 

1 = / ^^(27t)^6(k-pA+FA'), l = y -^^W^W-Pb+Pb'), 

1 = /^^^(27t)^5(ki+/i-;7i), \ = J-^^{2n)'^d{k3+h-pi): (A8) 
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/ dPi ! dp 1 dPi 



FIG. 9: The double gluon distribution. 


one can write the cross-section (IA6I) in a factorized form 

As^ /■ (f'k\(f'k2(f'k'} (f'hcf'hd^h 


da 


pp->ViV2X 


(27t)20 / k\klkl{ki +k2-kpf ijllllih +l2-h)^ 

X I J|3,<'™(P,,A:i,A:2,A:3,A:i+A:2-A:3) J r^a/$“|™(a/,/i,/2,/3,/i+/2 -/s) 

X / JaitJ|3i/S"*"Maifr,idi/,/ci i i ,z?i " (ait,Pi7,/c^ i ,Gi ,fi 


-5(A:i +/i —k^—h 


jdaipd^uS^^^^^aik, Pi/,^i±,^i±,Fi)5'i^^^*(oci;t5PiZ5^3±5^3±5j 

X J da2kd^2lS2^'^^{a2k,^2l,k2±,l2±,P2)S^2^'^*{a2k,^2l,k\ +k2-k3,li -f/2-/3,P2) 

X (27t)45(7ti -t- /i - Pi)d^p, {2%)^dik2 + l2-P2)d^p2 , (A9) 

and k 2 =k—ki,l 2 =l — h- The 2P -p impact factors from the third line of equation (IA9I) can be decom¬ 
posed into three color structures each, (see Fig.[9l), 

I = 5v:[8“'""8®“V(ri)s(ri -k2)f{k3)s{k3-kt) 

+ 8“'“38“J“V(*i)8(*i ~k2)f{k2)S(k2~ki) +8“'“-8“M/(ri)8(ri -r4)/(fe)8(fe-* 3 )]. 

each, and nine color terms are obtained after their multiplication. However, only one of these terms 
gives the leading contribution, the other ones being suppressed by the color factor (1 /N]:) or the limiting 
phase space {p\ ~ p 2 ). The normalization constant is determined comparing the collinear limit of 
equation (IA9I) to (IA2b . Indeed, using decomposition (lAlOl) and noting that the last three lines of (IA9I) 
convert into partonic cross-sections (IA4h . the differential cross-section (IA9I) can be written in the form 


d^app^v,V2X _ s^\e{N^ - 1)4 9(} 


da\d^\da2d^2 
dH 


(271) 


20 


X 


Then it follows that 


0ClPia2P2 

d^k 
(27t)- 


d^k\ 


f{auk\ 


7 


d^ki 

14 


f{a2,kl 


dH\ 

T 


/(Pi7i 


/ ^/(p2,/f)6i(ai,pi)a2(a2,p2) / 

(27t)2 


(AlO) 




.(A2-l)2’ 
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where the general eorrespondenee between the unintegrated f{a,k^) and eollinear g{(l,iS) gluon distri¬ 
butions 



was used. 
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